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ABSTRACT 

We report on our Chandra and RXTE observations of the bright old nova, 
V603 Aql, performed in 2001 April, supplemented by our analysis of archival 
X-ray data on this object. We find that the RXTE data are contaminated by the 
Galactic Ridge X-ray emission. After accounting for this effect, we find a high 
level of aperiodic variability in the RXTE data, at a level consistent with the 
uncontaminated Chandra data. The Chandra HETG spectrum clearly originates 
in a multi-temperature plasma. We constrain the possible emission measure 
distribution of the plasma through a combination of global and local fits. The 
X-ray luminosity and the spectral shape of V603 Aql resemble those of SS Cyg 
in transition between quiescence and outburst. The fact that the X-ray flux 
variability is only weakly energy dependent can be interpreted by supposing that 
the variability is due to changes in the maximum temperature of the plasma. 
The plasma density is likely to be high, and the emission region is likely to be 
compact. Finally, the apparent overabundance of Ne is consistent with V603 Aql 
being a young system. 

Subject headings: Stars: novae, cataclysmic variables — stars: individual (V603 
Aql) — X-rays: binaries 
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1. Introduction 

Cataclysmic variables (CVs) are interacting binary systems in which a white dwarf 
primary accretes from a Roche-lobe filling, late-type secondary (see Warner 1995 for a 
comprehensive review). In the majority of CVs, the magnetic field of the primary is not 
strong enough to control the accretion flow. In such systems (non-magnetic CVs), accretion 
proceeds via a Keplerian disk, which radiates half the available potential energy mostly in 
the optical and UV. The other half of the potential energy is in principle available to power 
the boundary layer between the disk and the white dwarf (Patterson & Raymond 1985a, b). 
It is the likely source of X-rays that are observed in nearby non-magnetic CVs, since the 
accretion disk is not hot enough to be a significant source of X-rays. There are, however, 
important questions regarding the validity of this simple interpretation. X-ray observations 
are an essential tool that may help advance our understanding of the accretion onto the 
white dwarfs in CVs in general, and the nature of the boundary layer in particular. 

A subset of CVs have been discovered as classical novae. These eruptions are due 
to thermonuclear runaway of hydrogen-rich material accreted on to the primary of a CV. 
V603 Aql was the brightest classical nova of the 20th century during its eruption (as Nova 
Aquilae 1918). Classical novae are spectacular laboratories of thermonuclear reactions, and 
make some notable contributions to the chemical evolution of the interstellar medium. How- 
ever, the outburst characteristics, such as the energetics, the duration, and nucleosynthesis 
yields, all depend on the details of the accretion processes in quiescence, including their 
quiescent mass accretion rate (M). Quiescent observations are therefore essential for the 
understanding of classical nova eruptions. 

Various authors have claimed V603 Aql to be magnetic; however, none of these claims 
have been verified by later, more sensitive, observations. The claim of periodicity in polar- 
ization by Haefner & Metz (1985) has been disproved by more sensitive observations by 
Cropper (1986) and by Naylor et al. (1996). Although there is a real modulation of the 
X-ray flux with a 60 minute timescale in in the Einstein data (Udalski & Schwarzenberg- 
Czerny 1989), this modulation is neither persistent or periodic, hence V603 Aql should be 
considered a non-magnetic CV (Borczyk et al. 2003; §4.1). 

V603 Aql is the brightest old (=quiescent) nova in the optical, UV and X-ray wave- 
lengths. It is also relatively nearby; Hubble & Duncan (1927) have estimated a distance 
d of 360 pc using the expansion parallax method. The statistical error is negligible in this 
estimate, although there is a systematic error due to departure from spherical symmetry in 
nova ejecta (Wade et al. 2000). Harrison & Gerz (1988) estimate d — 430 pc using the 
correlation between the absolute magnitude and the rate of decline of classical novae (t$ = 
8 days; Diirbeck 1988). At such a distance, Hipparcos parallax measurement does not lead 
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to a more precise estimate (d = 237^0° pc; Diirbeck 1999). In this paper, we adopt the 
Hubble & Duncan (1927) value of 360 pc as the fiducial distance. 

V603 Aql is also a permanent superhumper showing both positive and negative super- 
humps (Patterson et al 1997; Patterson 1999), suggesting a well-established disk with high 
M simultaneously undergoing apsidal and nodal precessions. Ferland et al. (1982) inferred 
a total disk luminosity of ~ 10 35 ergss -1 from UV data, requiring M > 10 18 gs -1 . V603 Aql 
has been observed in X-rays with Einstein , EXOSAT ROSAT PSPC, ASCA, and RXTE. 
Drechsel et al. (1983) found the source to be 3 times brighter in Einstein IPC in their ob- 
servation (~ 0.7 cts -1 , 1981 March) than in 1979 Oct (~ 0.3 cts -1 ). Drechsel et al. (1983) 
also inferred a very hard (kT>20 keV) spectrum from IPC and MPC data, and a total flux 
of L x = 3 x 10 33 egss -1 , a conclusion that we now consider suspect (§3.1). Note that even 
this luminosity can be generated by an accretion rate onto the white dwarf (M) of ~ 10 2 * * * * * * * * * * * * * 16 
gs -1 , much lower than that inferred from the accretion disk data. 

V603 Aql is the second brightest non-magnetic CV observed with ASCA (Baskill et 
al. 2004), which makes it a prime candidate for detailed X-ray spectroscopy. In this paper, 
we present the result of our Chandra X-ray Observatory High Energy Transmission Grating 
(HETG) observation of V603 Aql, with contemporaneous monitoring campaign using RXTE. 
We also briefly describe archival ASCA and RXTE data on V603 Aql. 


2. Observations 

V603 Aql was observed with ASCA (Tanaka et al. 1994) during 1996 October 8-10 

over a period of approximately 34 hrs. After standard screening, we have extaracted 64 s bin 

light curves in 2 bands (<2 keV and >2 keV) as well as over the entire ASCA energy range 

(0.5-10 keV), combining data from all 4 instrument. We have also extracted 2 spectra, one 

for the 2 GISs combined and the other for the 2 SISs combined, and associated response files. 

This observation is also featured by Baskill et al. (2004), who have analyzed the complete 

set of ASCA observations of non-magnetic CVs. 

V603 Aql was observed 3 times with RXTE (Bradt et al. 1993) during 1998 December 

5-7, each over a period of 6-7 hours, and over ~52 hours of total elapsed time. We have 

analyzed only the data obtained with the Proportional Counter Array (PCA; Jahoda et al. 

1996), which consists of 5 Proportional Counter Units (PCUs). We have applied the standard 

screening criteria, after which we have accumulated 41 ksec of exposure. The background 

was estimated using the faint source model using calibration file dated 2003 March 30. 

We have observed V603 Aql with Chandra (Weisskopf et al. 2002) using the High 
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Energy Transmission Gratings (HETG; Markert et al. 1994) and with the ACIS-S CCD 
array as detector (Garmire et al. 2003) on 2001 April 19/20 over a span of 18 hrs. Because 
Chandra is in a high orbit, there were no interruptions during the observation, resulting in 
an exposure of 63.5 ksec. We have used the pipeline-produced evt2 and pha2 files and used 
CIAO 3.0.2 with caldb 2.26 to generate the response file, which includes a calibration of the 
contaminant on the AGIS detector. We have combined ±lst orders of High Energy Grating 
(HEG) data into one spectral file, with the response file appropriately combined; similarly 
for Medium Energy Grating (MEG) data, with the ±lst orders combined. This observation 
was previously included in a comparative study of Chandra HETG spectra of 7 CVs (Mukai 
et al. 2003). 

Near the time of the Chandra observation, we have monitored V603 Aql with RXTE. 
We have obtained twenty short observations over 2001 April 17 through 22, each resulting 
in a typical exposure time of 2-3 ksec, for a total of 55.4 ksec of good time. For these as 
well as for the 1998 RXTE observations, light curves are presented in the form of ctss -1 per 
PCU, since variable numbers of PCUs were active at any given time. In addition, for the 
2001 observations, we have excluded the data taken with PCU0, as the propane leak in this 
detector makes background subtraction problematical 1 . 

We present summaries of these observations in Table 1. Earlier observations of V603 Aql 
with Einstein , EXOSAT, and ROSAT are summarized in Orio et al. (2001). 


3. Analysis and Results 
3.1. Contamination in Non-imaging X-ray Data 

V603 Aql is located in a crowded region near the Galactic plane (at l = 33.16°, b = 
0.83°). Background subtractions for RXTE PCA data are performed using a model that 
accounts for the time- variable particle background and for the average, high-latitude, cosmic 
X-ray background. Galactic diffuse emission remains in the “background subtracted” data, 
as well as any specific X-ray sources that are in the field-of-view. The latter includes the 
supernova remnant Kes 79 (Seward et al. 2003), 53.8' away, consisting of a thermal shell 
and a soft, central point source. This source is probably too soft and too distant to have a 
significant impact. More worrisome is the recurrent transient X-ray pulsar, GS 1843+009 
(Israel et al. (2001) and references therein), located 52.2' away from V603 Aql. During an 
outburst, this pulsar can be much brighter than V603 Aql, particularly at high energies since 


^ee http://rxte.gsfc.nasa.gov/docs/xte/oldnews/pcu-0.html. 
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GS 1843+009 is a hard source, and can strongly contribute to the counts of non-imaging 
observations aimed at V603 Aql. 


It is quite possible that the Einstein MPC data, taken in 1981 March and is reported by 
Drechse'l et ai. (1983), were strongly contaminated by G8 1843+009, which was unknown 
at the time (its discovery was in 1988; Koyama et al. 1990). Note that the MPC has a 
relatively wide (1.5° FWHM) filed-of-view and there are no contemporary all-sky monitor 
data. The high apparent ratio of MPC to IPC count rates (3.24/0.71) and the hard spectral 
shape inferred from the MPC data are both incompatible with our findings below. Either 
V603 Aql was in a different state in 1981 March, or the MPC data were contaminated by 
a hard source, such as GS 1843+009. We have extracted the MPC data from HEASARC 
and searched for the 29.5 s spin period of GS 1843+009; although we do not see the period, 
this unfortunately still allows the possibility of a significant contamination by this pulsar, 
since the pulse fraction of GS 1843+009 is unusually small (4% peak-to-peak; Koyama et al. 
1990). 


RXTE PCA is less susceptible to contamination by this source, because of its narrower 
field-of-view (1° FWHM). Moreover, we have checked RXTE ASM data to ensure that, at 
the time of these RXTE observations, GS 1843+009 was not in an active state. Although we 
cannot exclude contamination at a low level by GS 1843+009 or another source, it probably 
is not a major concern. 


However, the PCA data are contaminated by the apparently diffuse, “Galactic Ridge” 
X-ray emission (Warwick et al. 1985; Valinia & Marshall 1998). From Figure 1 of Warwick 
et al. (1985), we estimate a 2-6 keV flux of the Ridge of order 7xl0“ 12 ergss -1 cm -2 per 
EXOSAT ME beam (45' FWHM) near V603 Aql. Taking into account the larger beam size 
of RXTE PCA (1° FWHM), we expect of order 2 ctss -1 per PCU from the Ridge in the 
observations of V603 Aql 2 . Although the average spectrum of the Ridge has been precisely 
characterized, the fact that there are at least two components whose relative importance is 
a function of Galactic coordinates (Valinia & Marshall 1998) makes it difficult to subtract 
the Ridge contamination with confidence. In the rest of this work, we will use the RXTE 
data primarily to discuss the short timescale, relative variability of V603 Aql. 


2 Similarly, the Ridge definitely contributed significantly to the observed Einstein MPC count rate. How- 
ever, the Ridge spectrum is quite similar to that of V603 Aql, so a contribution from GS 1843+009 remains 
a likely possibility. 
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3.2. Variability 

We present the 64 s bin ASCA light (0.5-10 keV) in Figure 1. A significant variability 
is obvious: although the average count rate is 2.22 ct s” 1 , values ranging from <0.1 to ~7.0 
cts -1 (64 s bins, all 4 instruments and all photon energies combined) are seen. Despite the 
strong variability, there is no single stable period in this data (Baskill et al. 2004). We 
have also investigated the energy dependence of this variability. In Figure 2, we plot the 
softness ratio (counts below 2 keV divided the 2-10 keV counts) against the 2-10 keV count 
rate. Except at the lowest count rates, this plot is relatively flat, and shows a slight softening 
trend when the count rate is lower. This rules out changes in photoelectric absorption, which 
would predict a strong hardening at lower count rates, as a major cause of variability. 

We present the light curve of the 1998 RXTE observations in Figure 3. As the light 
curve contains counts due to the Galactic Ridge, we have chosen 2.0 (the estimated Ridge 
contribution) as the origin of the Y axis. With this offset, the high level of variability is 
apparent, including a large flare-like event at the beginning of the first observation. 

Although not the main strength, the Chandra HETG/ACIS-S observation provides an 
uninterrupted light curve of V603 Aql. Because the zeroth order (undispersed) image is 
heavily affected by pile-up, we restrict our analysis to the first order photons. We have 
combined the HETG 1st order light curve and the 2001 PCA light curves into a single figure 
(Figure 4). Of the 20 separate RXTE pointings, 4 overlapped in time with the Chandra 
observation. The middle panel shows the light curves of during this interval, both using the 
HETG (blue) and the PCA (black), the latter plotted with the ~2 cts -1 per PCU offset as 
before. These two light curves overlap closely. This suggests that our estimate of the Galactic 
Ridge contribution to the PCA data is a reasonable one. Moreover, any energy dependence 
of the variability must be small, since the two instruments are sensitive to different energy 
ranges. 

To investigate the energy dependence of the variability further, we have extracted spec- 
tra from the Chandra HETG (MEG) data from the high and low states, defined as when the 
total 1st order count rate was higher than/lower than 0.45 cts -1 (as shown by the dashed 
line in the middle panel of Figure 4). This ratio is shown as a function of the wavelength 
in Figure 5 in the 1.5-17.2A (0.72-8 keV) range. We confirm the ASCA result (Figure 2) of 
mild energy dependence. With this definition of high and low states, the high state spectrum 
is somewhat harder: a linear fit to the trend suggests high/low ratio is ~20% higher at 8 keV 
than at 1 keV. We reach a similar conclusion from the energy independence of the variability 
in the RXTE PCA data (without subtracting the Galactic Ridge, whose precise spectrum at 
the location of V603 Aql is uncertain). This suggests again that the variability of V603 Aql 
is only mildly energy dependent, and also that the spectra of V603 Aql and of the Galactic 
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Ridge are similar, barring an unfortunate coincidence. 

We present the power spectra (using the definition of Scargle 1982) from the 1998 
and 2001 observations in Figure 6. Baskill et al. (2004) contains a corresponding power 
spectrum plot from the ASCA observation. We confirm their finding that there is no coherent 
periodicity in the X-ray data of V603 Aql. In particular, the power spectrum of the Chandra 
data does not suffer from the aliasing problem that plagues low Earth orbit satellite data, 
and thus the multiple peaks seen in the bottom panel between 0.2-0. 5 mHz (periods in the 
2000-5000 s range) must be taken at face value. That is, V603 Aql shows strong variability 
on these timescales, which however is not due to a single, coherent period. 


In addition, we have attempted to verify the tentative detection of power at 60 hrs (Orio 
et al. 2001) using the 2001 RXTE campaign, but failed to do so. This is in part due to 
the the combination of the powers at higher frequencies (0. 2-0.5 mHz) and our sampling 
pattern. We also note that the earlier claim was based on ROSAT data, which leaves open 
the possibility of a soft X-ray absorption event that would be undetectable with RXTE. 
Finally, the highest peak in the Chandra power spectrum corresponds to a period of 25,000 
s. This is a significant fraction of the duration of our Chandra data, so we caution against 
overinterpreting this peak. 


3.3. Spectroscopy 

We show the average Chandra HETG/ACIS spectrum of V603 Aql in Figure 7. The line- 
rich nature of its spectrum is obvious. The He-like and cold Fe Ka lines are strongly detected, 
as well as a weaker H-like Fe Ka. The H-like lines of Si, Mg, Ne and O are prominent, while 
their He-like counterparts are detected but significantly weaker. In addition, many Fe L-shell 
lines are detected in the 10-17 A range. This overall appearance is that of an X-ray emission 
from multi-temperature, collisionally excited plasma, as was noted earlier by Mukai et al. 
(2003). 

We have analyzed this spectrum using two complementary methods. In one, we have 
characterized individual lines through local fits, using an unabsorbed power law continuum 
and Gaussian lines. In the second method, we have performed a global fit using the variable 
abundance version of the cooling flow model (vmkcf 1) with the Gaussian broadening function 
(gsmooth), as well as an additional Gaussian representing the fluorescent Fe line at 6.4 keV. 

In principle, we can derive several quantities from the global fit. Chief among them are 
the maximum temperature of the plasma, the global mass accretion rate, and the various 
elemental abundances. We report our global fit results in Table 2, with formal 90% confidence 
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errors. We do not detect a significant interstellar absorption; in fact, our upper limit (1.3 x 
10 19 cm -2 ) is considerably smaller than the value of 6-8xl0 20 cm -2 inferred from ROSAT 
PSPC data by Orio et al. (2001). This may indicate a possible problem with the calibration 
of contaminant for Chandra HETG/ACIS-S data, which would be a source of systematic 
errors in the low energy results, both from global and local fits. We also obtain only an 
upper limit for the minimum temperature, kT min . The fact that the model is pegged at its 
hard lower limit (0.0808 keV) suggests the limitation in the grid of models, in addition to 
the lack of sensitivity of HETG to such low temperature components, may be playing a role. 
As for the abundance values, we only report our results for Ne and Fe, two elements that 
are clearly non-solar relative to the fiducial values of Anders & Grevesse (1989). Other 
abundances are consistent with solar, although the best-fit values may deviate from 1.0 by 
0.1 or more. We consider the reliability and the implications of the global fit results further 
in §4.2 below. 

We list all the lines that are confidently detected in Table 3. For each line, we list 
the probable identification, observed wavelength, the velocity width, and the photon flux 
according to the local fit. These lines are also labelled in Figure?. The measured flux of 
a line can be combined with the theoretical curve of its emissivity in a collisionally excited 
plasma as a function of plasma temperature to draw a curve on the temperature - emission 
measure plane (Figure 8). For Ne and Fe, we have used the best-fit abundance values from the 
global fit in placing these curves. Any point on one of these curves represents the maximum 
emission measure a plasma of that temperature can have, before predicting more flux in that 
line than is observed. For multi- temperature plasma, emission measures must be summed 
at several discrete temperatures or integrated over the relevant range. This is analogous 
to the differential emission measure (DEM) distribution analysis favored for the study of 
stellar coronae (see, e.g., Griffiths & Jordan 1998). Note, however, we have not derived the 
true DEM distribution unambiguously. What we have done here is to constrain the shape 
of the DEM distribution from local fits, and to fit a discrete approximation of the DEM 
distribution appropriate for a isobaric cooling flow. Note also that our data for V603 Aql 
mixes lines from different elements, and lines of different types, and this complicates the 
interpretation. Nevertheless, we find it satisfactory that the line from the global fit parallels 
the locus defined by the series of parabola-like curves derived from individual line fluxes, 
with an appropriate offset (a factor of 4, which is roughly the number of temperature steps 
of the global model within the typical width of the parabola-like curves). 

It is also possible to study selected lines in greater detail. In Figure 9, we present a 
close-up view of the Fe Ka region of the spectrum as observed with the HEG. The equiva- 
lent width of the fluorescent (6.4 keV, or 1.98A) line is 150 eV, which is what is expected of 
reflection off a surface that subtends a 27r steradian as seen from the primary X-ray source, 
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if the Fe abundance is roughly solar. In Figure 10, we present close-up views of the He-like 
triplet regions for Mg and Ne. In both cases, we detect lines near the expected wavelength 
of the forbidden component. However, there appears to be an offset. Given this, the pos- 
sibility that these apparent detections are due to interlopers should be kept in mind. We 
place an upper limit to the flux of the forbidden component at laboratory wavelength of 
2.8xl0 -6 photons cm -2 s _1 for MgXI and 3.0xl0~ 4 * 6 photons cm~ 2 * * s _1 for NeIX, respectively. 
We explore the possible implications further in §4.4 below. 

The line width is generally small. Most lines are consistent with a velocity width of 
order 200-300 kms -1 . This is not surprising given the low inclination angle of V603 Aql 
(i ~ 20°; Patterson et al 1993). The exceptions are the Fe fluorescent line and several 
observed lines that are likely to be blends. 

Finally, we confirm that the same spectral model can be used to fit the ASCA spectra 
(those taken with GIS and with SIS) of V603 Aql. We obtain an acceptable fit just by 
adjusting the Nh and the participating mass accretion rate. The former is likely due to 
calibration issues (both for ASCA and for Chandra ), while the latter is not surprising, given 
the high level of variability known in this system. We obtain a participating mass accretion 
rate of ~ 1.6 x 10~ 10 M© yr -1 during the ASCA observation. We can also fit the Ridge- 
contaminated RXTE data with the same model, again suggesting the close similarity between 
the spectra of V603 Aql and the Ridge, but we refrain from further interpretation due to 
this problem. 


4. Discussion 

4.1. Confirmation of the non-magnetic nature of V603 Aql 

Our analysis of X-ray curves fails to show a persistent periodicity. Although the power 

spectra of individual observations often show prominent peaks, there is no single period that 

appears in all observations. Instead, the power spectrum of our Chandra data show the 

simultaneous presence of many peaks. Moreover, the nature of X-ray variability in V603 Aql 
is atypical for an IP. If there is an energy dependence in the X-ray variability of V603 Aql, it 
is rather subtle; a strong energy dependence, on the other hand, is a common characteristic 
of IPs. 

Combined with other recent X-ray results (Borczyk et al. 2003; Baskill et al. 2004), as 
well as non-detection of polarization (Cropper 1986; Naylor et al. 1996), we conclude that 
V603 Aql is a non-magnetic CV. The earlier results are likely correct in finding power at 
certain frequencies in the X-ray data, but are incorrect in interpreting the power spectrum 
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as evidence of a single, coherent, underlying clock. We consider the nature of the hard X-ray 
variability further in §4.3 below. 


4.2. Interpreting the Global Fit Results 

Scientifically, the fact that we must assume a physical model for the X-ray emission in a 
global fitting approach is both its strength and its weakness. It allows physical interpretations 
of the fit results, on the one hand. On the other hand, the interpretations are only as good as 
the set of assumptions one has adopted. In our global fits, we have assumed a isobaric cooling 
flow with no additional heating or cooling mechanisms, absorbed only by the interstellar 
medium. Several authors have already explored alternative assumptions while analyzing the 
Chandra HETG data of other CVs (Perna et al. 2003; Homer et al. 2004). Moreover, 
we cannot rule out the more simple model in which the emission measure is a power law 
function of temperature (cevmkl model in xspec; see, for example, Baskill et al. 2004). As 
Figure 8 shows, the isobaric cooling flow model deviates only mildly from a power law above 
temperature of 3xl0 6 K, where the HETG data have good sensitivity. In this respect, we note 
that our best-fit global model fails to reproduce the observed FeXVII line fluxes (Figure 7), 
and that a cevmkl model with a reduced O abundance is a viable model. Nevertheless, we 
discuss below our global fits using the isobaric cooling flow model, since it appears adequately 
to describe our data and has a physical interpretation. 

The numerical complexity and the degeneracy of model parameters are further factors 
limiting what one can achieve with global fits. For example, the maximum temperature 
(kT max ) can only be constrained if the overall abundance is fixed at an assumed value. This 
is because, for the likely range of kT max (>10 keV), the hottest plasma contributes little line 
fluxes, and the continuum shape within the Chandra HETG band changes only moderately. 
The main effect of increasing kT max is to reduce the equivalent widths of the lines, which 
can also result from lower abundances. We have therefore fixed the abundances to the solar 
values of Anders & Grevesse (1989), with the exception of Ne and Fe. Moreover, we have 
not included reflection in our fits, since Chandra HETG has little sensitivity to reflected 
continuum, even though the strength of the 6.4 keV line indicates that reflection must be 
occurring in V603 Aql. Given these assumptions, kT max is tightly constrained. Unlike in 
magnetic CVs, however, there is no obvious and simple relationship linking kT max with the 
white dwarf mass; this would require a knowledge of exactly how the plasma is shocked in 
an accretion disk boundary layer. 

The total mass accretion rate that participate in the cooling flow is well constrained, 
for the same solar abundance assumption, and for an assumed distance. However, as we 
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discuss in §4.4, some of the lines have the potential as plasma diagnostics; the other side 
of the coin is that the cooling flow model fit, which does not account for effects of line 
opacity, photoexcitation etc., may not be completely reliable. Moreover, the lines may have 
variable widths and/or profiles: physically, we expect the plasma to continue to decelerate 
as it cools, so the line width should be a function of the typical temperature needed to emit 
that line. We approximate this using the xspec convolution model, gsmooth, which allows 
the line width to be a function of the photon energy of that line. These obviously are not 
the same thing. However, the unconfused lines are marginally resolved at best in our data 
(see Table 3), so this is probably not a major concern. 

With these caveats, the mass accretion rate in the cooling flow is tightly constrained in 
our global fit (8.0 x 10 -11 M 0 yr -1 or 5.1 x 10 15 gs -1 ). In addition, the observed luminosity 
in the 1. 0-7.0 keV band is 1.88xl0 32 ergss -1 , depending only on the assumed d, and hence is 
more robust. That is, V603 Aql is hard X-ray luminous for a non-magnetic CV. In particular, 
V603 Aql is significantly more luminous than the dwarf nova Z Cam in outburst (6xl0 30 
ergss -1 ; Baskill et al. 2001) or the inferred hard X-ray component of the nova-like system, 
UX UMa (1.26xl0 31 ergss -1 ; Pratt et al. 2004). 

SS Cyg has a 3-20 keV luminosity of 3-9 xlO 31 ergss -1 during the peak of normal 
outbursts (Wheatley et al. 2003; McGowan et al. 2004), again fainter than V603 Aql. SS 
Cyg, however, has been observed at X-ray luminosity comparable to that of V603 Aql: in the 
1993 May ASCA observation obtained towards the end of an anomalous outburst (Nousek 
et al. 1994; Baskill et al. 2004) and transitions into and out of outburst (Wheatley et 

al. 2003; McGowan et al. 2004). Thus, SS Cyg has three main X-ray/optical states. In 
the presumed order of increasing mass accretion rate onto the white dwarf, these are: the 
quiescence, the transition, and the outburst peak. Of these three, V603 Aql resembles the 
transition most closely, in terms of its X-ray luminosity and its X-ray spectrum. The lack of 
detection of soft X-ray/EUV component in V603 Aql (Orio et al. 2001) is also consistent 
with this analogy. Arguing against this is the high optical/UV luminosity (10 35 ergs -1 ; 
Ferland et al. 1982) of its accretion disk. 

If the UV luminosity of V603 Aql indeed represents the accretion luminosity of its disk, 
the resemblance to the transition in SS Cyg must be a coincidence, and over 99% of the 
luminosity of the boundary layer remains unaccounted for. If this luminosity is radiated 
from the white dwarf surface, this must cover a large fraction of the white dwarf surface to 
have escaped detection: 10 35 ergs -1 radiated from 5% of the area of a 6,000 km radius white 
dwarf would have an effective temperature kT of 25 eV, and would have been detectable as 
a bright, soft component, even with a column of Nh = 8 x 10 2 ° cm -2 . At kT =20 eV, this is 
not necessarily the case, but the area must be 15% of the white dwarf surface. 
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An alternative interpretation is that the UV luminosity does not in fact arise from 
accretion, but instead due to irradiation by a hot white dwarf. Our understanding of when 
and how nuclear burning terminates in classical novae is still rather vague (see, e.g., Orio et 
al. (2001)), and little is known about the post-nova cooling. If the entire photosphere of 
the white dwarf in V603 Aql is hot, though less so than in super-soft source stage of classical 
novae, its photospheric emission can be hidden more easily than a soft X-ray component 
from a boundary layer, which must be hotter for a given total luminosity. In this picture, 
the hard X-ray luminosity represents a significant fraction of the accretion luminosity, and 
resemblance to SS Cyg in transition is not a coincidence. Further observations of the overall 
energy balance of SS Cyg during transition might provide a useful comparison to see if this 
model is feasible. 


4.3. The nature of the short time scale variability 

We have argued earlier that the rapid X-ray variability of V603 Aql is not due to 
photoelectric absorption. There is, however, one parameter in the cooling flow model that 
can explain the type of variability observed: kT max . When it is higher, the continuum is 
stronger and somewhat harder; when lower, weaker and somewhat softer 3 . In magnetic CVs, 
such changes would be unphysical, since the accretion flow is radial and kT max is strongly 
tied to the white dwarf mass. In a non-magnetic CV, however, the shock geometry is complex 
and three dimensional, and we know of no reason why the shock geometry, hence kT max , 
cannot vary with a timescale of 10s of minutes. On the other hand, this timescale is too long 
to be intrinsic to the boundary layer; some external mechanism to drive such a change is 
probably necessary. A search for similar X-ray variability in other non-magnetic CVs, such 
as SS Cyg during transition, would be worthwhile. 

This timescale, however, is similar to that seen in the UV spectroscopy of V603 Aql 
(Prinja et al. 2000). Based on the detection of episodic UV absorption lines, they suggested 
the presence of a time-variable wind in V603 Aql. This may be another symptom of the 
underlying mechanism that causes the X-ray variability. This is consistent with the fact that 
UV flux appears to be correlated with X-ray intensities (Borczyk et al. 2003). 


3 This does not contradict the earlier statement that the continuum shape does not strongly constrain 
kT max . The continuum shape can only constrain kT max to within ~10 keV, while the line to continuum 
ratio can constrain it to better than 1 keV, if the overall abundances are fixed. 


-13- 


4.4. Plasma Diagnostics 

The detection of individual line opens up the possibility to use selected line ratios as 
diagnostics of the temperature and the density of the emitting plasma. The first thing 
to note is that the high ratio of H-like to He-like line fluxes is a natural consequence of 
multi-temperature plasma, since H-like lines are typically emitted over a wider range of 
temperatures than are the He-like lines (see Figure 8). 

A casual glance at the He-like triplets of Mg and Ne (Figure 10) suggests that the 
forbidden component may be present, which would prove a relatively low density (n e < 
10 13 cm -3 ). However, we caution against this interpretation for several reasons. First, 
the observed wavelengths are not exactly at the laboratory wavelengths of the forbidden 
transitions, even though the resonant lines are exactly where they are supposed to be. This 
suggest either the identification is incorrect, or that the forbidden component originates in a 
material with different kinematic properties. Second, observed fluxes in the He-like triplets 
of Ne and O are much higher than those of the best-fit cooling flow model (Figure 7). This 
suggests the possibility of a second component. 

The third reason is that the strength of the fluorescent Fe Ka line (equivalent width ~150 
eV) is consistent with reflection from the white dwarf surface that subtends 27T steradians, 
as seen from the X-ray emitter. This strongly suggests that the emitting region is compact. 
The likely volume of such a compact boundary layer (much less than that of the white dwarf) 
and the emission measure estimates (Figure 8) require that the emission region must be of 
very high density. For a boundary layer covering 1% of the surface area of a white dwarf 
(with radius, r, of 6,000 km), and height O.lr, to have an emission measure of 10 54 cm -3 , 
the density must be 6 x 10 14 cm -3 . 

In fact, the lack of forbidden components of Mg XI and NeIX at rest wavelengths (see 
§3.3) is consistent with this interpretation. The ratio f/(r+i) is measured to be <0.44 for 
Mg XI and <0.16 for NeIX, taking the upper limit for the / flux and the lower limits for r 
and i. This implies density greater than ~ 4 x 10 13 cm -3 (Porquet & Dubau 2000), and can 
be much higher. 

What about the Fe L diagnostics (Mauche et al. 2001, 2003)? Of the FeXVII lines, the 
17.05A line appears to be strongly detected, while the 17.10A line is not confidently detect- 
ed. Within the limited S/N of our data, this is consistent with high density or significant 
photoexcitation. As for the FeXXII lines near 11. 8A the S/N of our data is not high enough 
to allow reliable use of this diagnostic, given the need to first separate the FeXXII 11. 78 A 
line from the Fe XXIII 11.74A line. 
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4.5. Abundances 

Our analysis suggests that the accreting gas in V603 Aql has a high Ne abundance and 
a low Fe abundance, relative to the fiducial values of Anders & Grevesse (1989). However, 
note that the solar photospheric Fe abundance more recently derived by Lodders (2003) is 
approximately 60% that of Anders & Grevesse (1989). Thus, whether the Fe abundance 
of V603 Aql is solar or sub-solar may well be a question of the Fe abundance in the Sun. 
However, the Ne overabundance of V603 Aql appears genuine, regardless of which solar 
values we adopt. 

V603 Aql appears to differ from other CVs for which abundance anomalies have been 
reported (see, for example, Mauche et al. 1997; Gansicke et al. 2003 and references therein), 
since these concern mostly carbon and nitrogen, and are likely caused by an episode of CNO 
cycle in these binaries. Although further nuclear evolution can in principle alter the Ne 
abundance, this presumably would have been preceded by the CNO cycle. Since the NV 
A1240 line is seen to be weak (Ferland et al. 1982), the CNO cycle and further nuclear 
evolution does not appear to be a viable explanation for the Ne overabundance in V603 Aql. 

Approximately a third of all classical novae are thought to harbor an O-Ne-Mg white 
dwarf, as evidenced by the high Ne abundances in their ejecta (Starrfield et al. 1986; 
Politano et al. 1995). The X-rays are emitted just as matter accretes onto the white dwarf 
surface, and hence likely represents the current composition of the secondary. If V603 Aql is 
a O-Ne-Mg nova, the pollution of the secondary by the Ne-rich ejecta is possible in principle. 
However, according to the analysis of Pottasch (1959), the Ne abundance of the V603 Aql 
shell was 1/7,000, or 1.4xl0~ 4 by number. This value is between the Solar value of Anders 
& Grevesse (1989) (1.23xl0~ 4 ) and our X-ray fitting result (1.75xl0~ 4 ). We conclude that 
V603 Aql is not a O-Ne-Mg nova. 

This suggests that the V603 Aql system was born with a higher then solar Ne abundance 
due to Galactic chemical evolution, implying that it is either young, or it originated closer to 
the Galactic center (see, e.g., Chiappini et al. 1997). Given the location of V603 Aql very 
near the Galactic mid-plane ( b = 0.83° at a distance of 360 pc), the former interpretation is 
certainly plausible. It is also consistent with the proposition that there are two population of 
classical novae, one young and fast (V603 Aql being a prominent example), the other older 
and slower (Della Valle et al. 1992). Abundance measurements of other old novae may 
provide a way to test whether there indeed are two distinct populations. 
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5. Conclusions 

We have presented our analysis of Chandra, RXTE, and archival X-ray observations 
of V603 Aql. It is among the X-ray brightest CV and displays a line-rich X-ray spectrum. 
However, it is likely that non-imaging X-ray data of V603 Aql are contaminated by the 
Galactic Ridge as well as nearby, unrelated sources. On V603 Aql itself, our findings can be 
summarized as follows. 

(1) It is a highly variable X-ray source. However, we find no indications of a single, 
coherent period, and confirm that V603 Aql is a non-magnetic CV. 

(2) The X-ray spectrum of V603 Aql is from a multi-temperature, collisionally ionized 
plasma. Specifically, we find that the cooling flow model provides an adequate framework 
for physical interpretation. 

(3) The spectrum and the luminosity of V603 Aql resembles those of SS Cyg in transition, 
even though the UV luminosity of V603 Aql is very high. One possible interpretation is that 
the accretion rate in V603 Aql is similar to that of SS Cyg in transition, even though this 
leaves the question of the origin of the high UV luminosity unanswered. 

(4) The weak energy dependence of the X-ray flux variations in V603 Aql can be ex- 
plained as due to changes in the maximum temperature of the multi-temperature plasma. 
The physical cause of such changes is unknown, but we note the accretion disk wind in 
V603 Aql appears to be modulated on a similar timescale. 

(5) There are complications that make the plasma diagnostics less clear-cut than we 
would like. The preponderance of evidence points towards a high density, compact emission 
region compatible with the accretion disk boundary layer. 

(6) We detect an apparent overabundance of Ne in the accreting matter. This is best 
explained as due to V603 Aql being a young binary. 


This research has made use of data obtained from the High Energy Astrophysics Science 
Archive Research Center (HEASARC), provided by NASA’s Goddard Space Flight Center. 
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Fig. 3. — RXTE PCA light curves of in 16 s bins from the 1998 observations. Particle and 
high Galactic latitude background have been subtracted, but not the Galactic background. 
The latter is estimated to be approximately 2 cts -1 , which has been chosen as the origin of 
Y axis of this figure. 
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Fig. 4. — RXTE PCA light curves in 16 s bins, and Chandra HETG 1st order light curves in 
60 s bins, of V603 Aql from the 2001 observations. The middle panel shows the HETG curve 
in blue (with scales labelled on the right) and the PCA curves from 4 overlapping pointings 
in black. The horizontal dashed line is the HETG count rate used to divide the data into 
high and low states for Figure 5. The top and bottom panels show the non-overlapping PCA 
curves; all PCA light curves are offset by 2 cts -1 . 
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Fig. 5. Ratio of the Chandra HETG (MEG) high state spectrum to the low state spectrum 
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Fig. 6.— Power spectra of RXTE PCA and Chandra HETG data plotted against frequency 
in mHz. 
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Fig. 7. — The average spectrum of V603 Aql as observed with Chandra HETG/ACIS-S, 
plotted in 3 sections. For clarity, we plot only the HEG data in the top panel, and only the 
MEG data in the bottom two panels. The data are shown as black histograms, while the 
best fit cooling flow model (see text for details) are shown as red lines. Emission lines that 
are significantly detected are labeled. 
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Fig. 8. — The emission measure distribution as derived from the global fit, and the constraints 
on it derived from individual line fits. The latter curves are shifted down by a factor of four, 
which approximately accounts for the width of these parabola-like curves compared with the 
sampling of the emission measure distribution from the global fit. Blue curves are for H-like 
lines, red for He-like lines; black dashed lines are for Fe L shell lines. 
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Table 1. Log of Observations 


Observatory 

Observation ID 

Date 

Exposure 

(ksec) 

Average Count 
(cts -1 ) 

ASCA 

34013000 

1996 October 8-10 

38(GIS)/ 43(SIS) 

2.22 (total) 

RXTE 

30025-01“ 

1998 December 5-7 

41 (total) 

3.00 b 

RXTE 

60013-01 c 

2001 April 17-22 

55.4 (total) 

3.25 b 

Chandra 

1901 

2001 April 19/20 

63.5 

0.33 d 


a Full Observation IDs are 30025-01-01-00, 30025-01-01-01, and 30025-01-01-02. 

b Per PCU. Includes contribution from the Galactic Ridge, estimated to be ~2 cts -1 per 
PCU. 

c Full Observation IDs are of the form 60013-01-nn-00, where nn is in the range 01 to 20. 
d ±lst orders combined. 


Table 2. Global Fit Results 


N h 

kTmin 

kT 

^ ± max 

Ne a 

Fe a 

M h 

< 1.3 x 10 19 cm -2 

< 0.14 keV 

20.31} i keV 

1421 9 ; 21 

0.621°$ 

8.00±0.01 x 10" 11 M, 


“Abundances, relative to solar values of Anders & Grevesse (1989) 

b Participating mass accretion rate, for a fiducial distance of 360 pc (Hubble & Duncan 
1927). 
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Table 3. Emission Lines 


Line ID 

Wavelength 

(A) 

Width 
(km s _1 ) 

Flux 

(10 -6 photons cm _2 s' ! ) 

Fe XXVI 1.78 

1.774 (1.767-1.787) 

1150 (<3050) 

12.5 (4.9-24.2) 

FeXXV 1.86 

1.855 (1.851-1.860) 

550 (<2250) 

24.3 (15.2-40.3) 

Fe Fluorescent 

1.939 (1.933-1.945) 

1520 (770-2680) 

20.9 (12.5-31.3) 

SXVI 4.73 

4.732 (4.713-4.738) 

160 (<660) 

5.9 (2.6-10.1) 

SXV 5.04 

5.043 (5.034-5.048) 

0 (<800) 

4.9 (2.8-10.5) 

Si XIV 6.18 

6.182 (6.180-6.184) 

300 (180-430) 

19.3 (16.4-22.7 ) 

Si XIII 6.65 

6.651 (6.647-6.659) 

460 (230-970) 

12.0 (8.7-15.7) 

Fe XXIV 7.99 

7.981 (7.976-7.990) 

370 (10-750) 

3.8 (2. 0-6.0) 

Mg XII 8.42 

8.420 (8.418-8.422) 

280 (200-400) 

18.0 (15.3-21.1) 

Mg XI 9.17 

9.170 (9.167-9.175) 

0 (< 180) 

5.7 (4.0-8.4) 

Mg XI 9.23 

9.230 (9.220-9.237) 

0 (< 180) 

2.8 (1.2-5. 2) 

Mg XI 9.31 ? a 

9.335 (9.326-9.340) 

0 (<600) 

3.2 (1. 1-5.9) 

FeXX 10.00 

9.978 (9.975-9.988) 

0 (<370) 

6.6 (4. 1-9.6) 

FeXX 10.05 

10.051 (10.041-10.061) 

0 (<360) 

2.8 (1. 1-5.7) 

NeX 10.24 

10.245 (10.240-10.254) 

0 (<340) 

6. .4 (3. 6-9. 7) 

Blend b 

10.644 (10.637-10.653) 

660 (480-990) 

24.0 (18.0-30.8) 

Blend c 

11.019 (11.010-11.029) 

590 (350-820) 

18.7 (13.4-25.1) 

FeXXIV 11.18 

11.177 (11.170-11.183) 

400 (230-690) 

16.1 (11.4-21.9) 

Fe XXIV 11.43 

11.439 (11.429-11.440) 

0 (<280) 

10.1 (6.8-15.3) 

Blend d 

11.756 (11.737-11.773) 

820 (510-1280) 

15.4 (8.3-23.8) 

NeX 12.13 

12.135 (12.132-12.139) 

270 (120-360) 

45.6 (36.2-53.4) 

NeIX 13.45 

13.439 (13.431-13.445) 

220 (100-380) 

19.8 (14.2-30.7) 

NeIX 13.55® 

13.531 (13.516-13.544) 

800 (580-1130) 

37.3 (27.2-48.3) 

NeIX 13.70 ? f 

13.744 (13.730-13.756) 

70 (<570) 

7.5 (2.5-12.3) 

FeXVII 14.37 

14.360 (14.355-14.368) 

0 (<280) 

8.8 (4.7-16.1) 

FeXVII 15.01 

15.009 (15.001-15.014) 

270 (10-480) 

46.3 (23.7-48.9) 

FeXVII 15.26 

15.260 (15.256-15.267) 

0 (<190) 

15.5 (11.3-26.5) 

Fe XVIII 15.63 

15.623 (15.607-15.641) 

300 (100-630) 

8.0 (4.4-20.7) 

Fe XVIII 16.05 

16.035 (15.996-16.054) 

760 (460-1040) 

31.3 (19.2-45.4) 

FeXVII 16.78 

16.780 (16.764-16.789) 

370 (260-650) 

30.6 (19.2-44.0) 


Table 3 — Continued 


Line ID 

Wavelength 

Width 

Flux 


(A) 

(km s -1 ) 

(10 -6 photons cm -2 s -1 ) 

FeXVII 17.05 

17.061 (17.052-17.071) 

440 (340-650) 

49.9 (35.5-67.1) 

OVIII 18.97 

18.964 (18.958-18.970) 

330 (270-430) 

144.2 (115.6-169.3) 


a Flux at the rest wavelength of the forbidden Mg xi is < 2.8 x 10 6 photons cm~ 2 s~ 1 

b Blend of FeXIX 10.64 and FeXIX 10.65? 

c Blend of Fe XXIII 11.02 and FeXXIV 11.03? 

d Blend of FeXXIII 11.74 and FeXXII 11.77? 

e Possibly blended with Fe,XIX 13.52 

f Flux at the rest wavelength of the forbidden Ne ix is < 3.0 x 10~ 6 photons cm -2 s -1 



